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Beschreibung 

Identif izieren pharmazeutischer Targets 

Die menschliche Erbsubstanz (Genom) umfasst schatzungsweise 
20.000 bis 80.000 Gene, die den genetischen Code fur etwa ei- 
ne Million EiweiJistof f e (Proteine) beinhalten. In den spezia- 
lisierten Korperzellen werden jeweils nur Untermengen aller 
Gene tatsachlich abgelesen (exprimiert ) . Die Gesamtheit der 
dadurch erzeugten Proteine wird als Proteom dieser Zelle be- 
zeichnet. Das Wechselspiel der Proteine untereinander sowie 
mit der DNA stellt den wichtigsten Teil der Maschinerie dar, 
die der Entwicklung des menschlichen Korpers aus der befruch- 
teten Eizelle sowie alien Korperf unktionen zugrunde liegt. 
Aus der Sicht der Informatik stellt die Erbsubstanz damit ei- 
nen prozeduralen Code fur die Struktur und Funktion des 
menschlichen Korpers dar. 

Viele Krankheiten und Fehlf unktionen des Korpers gehen auf 
Storungen des f unktionellen Netzwerks aus Genom und Proteom 
zurlick. So wirken einige Medikamente als Agonisten bzw. Anta- 
gonisten spezifischer Zielproteine, d. h. sie verstarken oder 
schwachen die Funktion eines Proteins mit dem Ziel, das aus 
Proteom und Genom gebildete regulatorische Netzwerk zuriick in 
einen normalen Funktionsmodus zu bringen. Diese Zielproteine 
(Targets) werden bislang nach heuristischen Prinzipien aus 
biochemischen Uberlegungen abgeleitet. Oft ist dabei unklar, 
ob die Fehlfunktion eines Proteins tatsachlich die Krank- 
heitsursache oder nur eines der Symptome einer versteckten 
Fehlregulation an anderer Stelle des Netzwerks darstellt. 

Aus [1] und [5] ist bekannt, Nervenzellen (Neuronen) und de- 
ren biologische Funktionalitat, d.h. deren biologisches Ver- 
halten, durch kunstliche Neuronen nachzubilden . 

Ferner ist aus [1] bekannt, zwei Arten von Synapsen bzw. Ner- 
venzellen zu unterschieden, eine erregende (exzitatorische) 
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Synapse 251 bzw. Nervenzelle und eine hemmende (inhibitori- 
sche) Synapse 252 bzw. Nervenzelle. 

Hemmende Synapsen 252 senken zu ubertragende bzw. weiterzu- 
5 leitende elektrische Potentiale, erregende Synapsen 251 erho- 
hen zu ubertragende elektrische Potentiale. 

Weitere Ausfuhrungen zum Aufbau und zur Funktionalitat einer 
Nervenzelle sowie zur Nervenleitung sind in [1] gegeben. 

10 

Ferner ist aus [1] und [5] eine kiinstliche Nervenzelle 
(kiinstliches Neuron) , welche eine (biologische) Nervenzelle 
nachbildet bekannt. 

15 Anschaulich gesehen ist ein solches kunstliches Neuron eine 
mathematische Abbildung, welche entsprechend dem Ubertra- 
gungsverhalten der biologischen Nervenzelle eine Eingangsgro- 
£ e d. es kiinstlichen Neurons auf eine Ausgangsgrofie des kiinst- 
lichen Neurons abbildet. 

20 

In Anlehnung an das biologische Vorbild besteht ein kunstli- 
ches Neuron aus drei Komponenten : dem Zellkorper, den Dendri- 
ten, welche Eingangssignale in das kiinstliche Neuron aufsum- 
mieren, dem Axon, welches das Ausgangssignal des kiinstlichen 
25 Neurons nach auJien weiterleitet, sich verzweigt und mit den 

Dendriten nachf olgender kiinstlicher Neuronen iiber Synapsen in 
Kontakt tritt. 

Eine Starke einer Synapse bzw. die Art der Synapse wird meist 
30 durch einen numerischen Wert bzw. dessen Vorzeichen darge- 
stellt. Dieser Wert wird als Verbindungsgewicht bezeichnet. 

Entsprechend dem biologischen Vorbild lasst sich ein Ubertra- 
gungsverhalten bzw. das Abbildungsverhalten eines kunstli- 
35 chen Neurons wie in [1] und [5] beschrieben abbilden. 
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Weitere Ausfiihrungen zu kiinstlichen Neuronen und deren Funk- 
tionalitat sind in [1] und [5] gegeben. 

Ferner ist aus [2] und [5] bekannt, einzelne Neuronen mitein- 
5 ander zu verknupfen. Ein solche Anordnung miteinander ver- 

kniipfter Neuronen wird als neuronales Netz bezeichnet. Grund- 
lagen iiber neuronale Netze, beispielsweise verschiedene Arten 
von neuronalen Netzen, Trainingsverf ahren fur neuronale Net- 
ze, Beziige zu biologischen Nervenzellen-Anordnungen, sind in 
10 [2] beschrieben. 

Aus [3] und [5] ist eine Anwendung eines Mean-Field Modells 
bei der Beschreibung eines komplexen Systems bekannt. Bei ei- 
nem Mean-Field Modell werden stochastische Wechselwirkungs- 
15 einflusse zwischen Komponenten eines Systems durch einen 

mittleren Wechselwirkungseinf luss genahert. Damit lassen sich 
analytisch nicht beschreibbare, stochastische Systeme auf be- 
schreibbare, deterministische Systeme reduzieren. 

20 Aus [4] und [5] ist die Anwendung des Mean-Field Modells bei 
der Beschreibung einer Neuronenstruktur bekannt. 

) Aufgabe der Erfindung ist es, das Identif izieren von Protei- 

nen, die sich als Ziel medikamentoser Behandlung genetisch 
25 bedingter Krankheiten oder Storungen eignen, zu verbessern. 

Diese Aufgabe wird durch die Erfindungen gemaii den unabhangi- 
gen Anspriichen gelost. Vorteilhafte Weiterbildungen der Er- 
findungen sind in den Unteransprilchen gekennzeichnet . 

30 

Erf indungsgemafi wird eine Mehrzahl von Genexpressionsmustern 
gleichartiger Zellen oder eines Gewebes bestimmt, wobei je~ 
weils die Expressionsrate der Gene der Zelle bestimmt wird. 
Die Mehrzahl der Genexpressionsmuster wird -derart bestimmt, 
35 dass der zeitliche Verlauf der Genexpressionsmuster der Zelle 
zumindest teilweise rekonstruiert werden kann. Ein dynami- 
sches Modell des regulatorischen Netzwerks aus Genom und Pro- 
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teom der Zelle wird gebildet, indem ein aquivalentes Neurona- 
les Netz in der folgenden Weise gebildet wird: 

i) Ein Gen des Genoms sowie das zugehorige Protein werden 
durch ein Neuron des aquivalenten Neuronalen Netzes reprasen- 

5 tiert. 

ii) Die Expressionsrate eines Gens wird durch eine nicht- 
negative Aktivitat des aquivalenten Neurons reprasentiert. 

iii) Die regulatorische Wirkung eines Proteins auf ein 
Gen wird durch eine synaptische Verbindung von dem zum Prote- 

10 in aquivalenten Neuron zu dem zum Gen aquivalenten Neuron 
reprasentiert . 

iv) Die Art einer regulatorischen Wirkung (verstarkend 
oder inhibitorisch) wird im Neuronalen Netz durch das Vorzei- 
chen und die Starke eines zugehorigen synaptischen Gewichts 

15 reprasentiert. 

v) In einer Weiterentwicklung ist es auch moglich, eine 
posttranslationale Modifikation eines ersten Proteins durch 
ein zweites Protein durch eine synaptische Verbindungen mit 
rnultiplikativer Wirkung vom zweiten Neuron zum ersten Neuron 

20 zu reprasentieren. 

vi) In einer weiteren Weiterentwicklung kann ein externer 
Einfluss auf das regulatorische Netzwerk durch einen Input- 
knoten des aquivalenten Neuronalen Netzes reprasentiert wer- 
den.. 

25 

Das aquivalente Neuronale Netz wird mit den bestimmten Ge- 
nexpressionsmustern verglichen und an diese angepasst. Aus 
dem angepassten Neuronalen Netz wird das regulatorische Netz- 
werk der untersuchten Zelle erschlossen. 

30 

Es wird also eine Aquivalenzbeziehung zwischen dem funktio- 
nellen Netzwerk des Genoms und Proteoms einerseits und dem 
Neuronalen Netz des menschlichen Gehirns andererseits, die 
beide stark vernetzte riickgekoppelte Systeme darstellen, her- 
35 gestellt. Durch diese Abbildung gelingt die Modellierung des 
funktionellen Netzwerks aus Proteinen und Genen. 
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Das Verfahren erlaubt die Identif izierung von Targetproteinen 
auf systematischer Basis. Zwischen genetischen und Neuronalen 
Netzen kann die beschriebene Aquivalenzbeziehung hergestellt 
werden. Das dynartiische Wechselspiel von Genen und regulatori- 
5 schen Proteinen wird damit durch ein dynamisches Neuronales 
Netz modelliert. Das Verfahren nutzt die in der Zeitabfolge 
der Genexpressionsmuster enthaltene Information fur die Iden- 
tif izierung kausaler regulatorischer Zusammenhange. 

10 In der Regel fuhren genetisch bedingte Krankheiten zu Tcomple- 
xen Fehlfunktionen, die aber oft auf nur ein oder wenige feh- 
lerhaft arbeitende Gene oder Proteine zuruckgehen . Bislang 
sind diese Schliisselgene bis auf Einzelfalle nicht bekannt . 
Stattdessen werden im herkommlichen Target Finding durch heu- 

15 ristisches Vorgehen Zielproteine (Targets) gesucht, deren me- 
dikamentose Regulierung das gesunde Genexpressionsmuster 
bestmoglich wiederherstellen soil. 

Neuartige Schatzungen sprechen von etwa 10.000 verschiedenen 
20 Proteinen als mogliche Targets im menschlichen Genom, deren 
Durchf orstung durch einen heuristischen Ansatz allein nicht 
zu bewerkstelligen ist. 

Der hier beschriebene Modellansatz stellt ein leistungsf ahi- 
25 ges Verfahren zum systematischen Target Finding dar, also zur 
Identif ikation eines oder weniger Schliisselgene oder 
-proteine, die am Anfang der Regulationskaskade stehen und 
die beispielsweise eine Organentwicklung einleiten, die Rege- 
nerationsf ahigkeit von Gewebe regulieren oder aber auch fur 
30 raeist komplexe Veranderungen von Genexpressionsmustern im 
Krankheitsf all verantwortlich zeichnen. 

Das beschriebene Verfahren erlaubt ein computerunterstiitztes 
Target Finding, das in der Lage ist, die grolie Datenmenge und 
35 die vielfaltigen und komplexen Zusammenhange zu analysieren. 
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Es lassen sich damit folgende Anwendungsgebiete spezifizie- 
ren : 

- Modellbasierte Unterstiitzung von Forschungsaktivitaten zur 
5 Entschliisselung der menschlichen Morphogenese und damit der 

generellen Prinzipien von genetisch gesteuerten Wachstums-, 
Regenerations- und Abbauprozessen. 

- Unterstiitzung der Identif ikation von Zielproteinen, die fur 
10 Wachstumsstorungen wie z. B. ungehemmtes Tumorwachstum ur- 

sachlich verantwortlich sind und nicht nur ein Symptom dar- 
stellen. Daraus konnen neuartige Verfahren fur eine hochsen- 
sitive Tumorf riihdiagnostik aber auch Behandlungsmethoden wie 
der selektiv induzierte Zelltod (Apoptose) in Tumorzellen ab- 
15 geleitet werden. 

- Unterstiitzung der Identif ikation von regulatorischen Prote- 
inen, die in Wachstums- und Regenerationsprozesse eingreifen. 
Damit ware eine wichtige Htirde auf dem Wege der induzierten 

20 Regeneration von Organen genonunen. 

Im Folgenden wird die Erfindung anhand von Ausf uhrungsbei- 
spielen naher erlautert, die in den Figuren schematisch dar- 
25 gestellt sind. Gleiche Bezugszif f ern in den einzelnen Figuren 
bezeichnen dabei gleiche Elemente . Im Einzelnen zeigt: 

Fig. 1 schematisch die regulatorischen Vorgange, die das 



30 



Fig. 2 



Expressionsmuster einer Zelle bestimmen; 

eine schematische Darstellung der Vernetzung von 



Fig. 4 



Fig. 3 



Neuronen; 

die Potenziale innerhalb eines Dendriten bzw. eines 
Neurons als Funktion der Zeit; und 

eine schematische Darstellung einer modulatorischen 



35 



Synapse . 
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Fig. 1 zeigt die wichtigsten Wechselwirkungen zwischen Genen 
und Proteinen eines DNA-Abschnitts auf . Die Wechselwirkungen 
werden als Basis fur die Beschreibung des genomischen regula- 
torischen Netzwerks herangezogen. 

5 

Im oberen Teil der Fig. 1 ist schematisch angedeutet, wie ein 
von aufien auf die Zelle einwirkendes externes Signal - etwa 
im Rahmen der interzellularen Kommunikation das beispiels- 
weise von einem Transmembran-Rezeptorprotein (z. B. von einem 
10 Kalziumkanal) aufgenommen und in geeigneter Weise in das In- 
nere der Zelle ubertragen wird, die Produktion der Gene A, B, 
C und D des DNA-Abschnitts auslost. 

Es besteht daher prinzipiell auch die Moglichkeit, die Ex- 
15 pressionsrate einzelner Gene einer Zelle uber die^ erwahnten 
Wege von aufierhalb der Zellen zu beeinf lussen . 

Als ein Gen wird ein nicht notwendigerweise zusammenhangender 
Abschnitt der DNA bezeichnet, der den genetischen Code fur 

20 ein Protein oder auch fur eine Gruppe von Proteinen enthalt . 
Im allgemeinen weist die DNA sogenannte Exons und Introns 
auf. Exons stellen Teile der DNA dar, die tatsachlich ein 
Protein codieren. Introns stellen Teile der DNA dar, die 
nicht unmittelbar ein Protein codieren. In erster Naherung 

25 sind sie ohne Funktion. Exons und Introns wechseln einander 
in der DNA ab. Bezeichnet man als ein Gen die Menge der E- 
xons, die zusammen ein bestimmtes Protein codieren, so ist 
ein solches Gen - wie oben erwahnt - in der Regel nicht zu- 
sammenhangend . 

30 

Der Produktionsvorgang eines Proteins aus einem Gen, zum Bei- 
spiel Protein A ausgehend von Gen A in Fig. 1, wird als Ex- 
pression dieses Gens bezeichnet. Die Ubersetzung des DNA- 
Codes des Gens in die Kette der Aminosauren des Proteins wird 
35 als Translation bezeichnet. Die Rate, mit der das Protein A 
in einem gegebenen Kontext produziert wird, wird seine Ex- 
pressionsrate genannt. 
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Nicht alle Gene werden in einer Zelle exprimiert . Vielmehr 
unterscheiden sich verschiedene Zelltypen durch ihr Ge- 
nexpressionsmuster . Dies gilt oftmals auch fur den Unter- 
schied zwischen kranken und gesunden Zellen. 

Das Expressionsmuster einer Zelle wird durch die in Fig. 1 
schematisch dargestellten regulatorischen Vorgange bestimmt. 
Die regulatorischen Vorgange werden im Wesentlichen von eini- 
gen wichtigen Wechselwirkungen zwischen Proteinen und Genen 
sowie zwischen den Proteinen untereinander bestimmt. 

So kann die Expressionsrate eines Gens A durch das Vorhanden- 
sein eines anderen Proteins B reguliert, d. h. erhoht, er- 
niedrigt oder zum Erliegen gebracht werden. In diesem Bei- 
spiel wirkt das Protein B regulatorisch auf das Gen A bzw. 
das Protein A. Zu regulatorischen Proteinen konnen beispiels- 
weise die Proteinbausteine von Aktivatorkomplexen gerechnet 
werden. Regulatorische Proteine konnen sich gleichzeitig auf 
viele Zielgene auswirken. 

Eine zweite Art der Wechselwirkung besteht in der posttrans- 
lationalen Modifikation von Proteinen, d. h. der Modifikation 
von Proteinen nach der Translation. In der Regel erfolgt die 
posttranslationale Modifikation eines Proteins im unmittelba- 
ren Anschluss an die Translation, d. h. bevor das Protein in 
der Zelle wirkt. So werden zum Beispiel viele Proteine von 
speziellen Enzymen phosphoryliert oder glykolysiert, d. h. 
das Zielprotein wird durch Anhangen bzw. Abspalten chemischer 
Gruppen in seinen f unktionellen Zustand gebracht oder in ei- 
nen Zustand versetzt, in dem es nicht mehr wirksam ist. 
Posttranslationale Modifikation kann also ein Protein gegebe- 
nenfalls temporar funktionell einschalten oder ausschalten. 

In Fig. 1 ist das Protein A ein sogenanntes Ef f ektorprotein, 
d. h. es wirkt innerhalb der Zelle auf andere Substanzen und 
nicht unmittelbar auf das Genom oder Proteom. In Fig. 1 modi- 
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fiziert das Protein C im Wege cler posttranslationalen Modifi- 
kation die Funktion des Ef f ektorproteins A. 

Protein B ist ein regulatorisches Protein, da es die Expres- 
5 sionsrate des Proteins A bestimmt, indem es mit demjenigen 

DNA-Abschnitt wechselwirkt, der das Gen A enthalt . Das Prote- 
in D modif iziert somit die Funktion eines regulatorischen 
Proteins (Protein B) im Wege der posttranslationalen Modifi- 
kation. 

10 

Datenbasis 



Die Nukleinsauresequenz der menschlichen DNA ist weitestge- 
hend bekannt. Auch die durch die DNA codierten Gene sind in 

15 zunehmendem Mafie identif iziert . Nicht ganz so vollstandig ist 
das Wissen iiber das Proteom, einschlieftlich der eventuell 
durch Wechselwirkung zwischen den Proteinen posttranslational 
modif izierten Proteine. Allerdings erlauben neuere Sequenzie- 
rungs- und Hochdurchsatz-Screeningverf ahren eine rasche Iden- 

20 tifizierung weiterer Gene und Proteine. 

Ein weiterer wichtiger Schritt zur Aufklarung der Expressi- 
onsmuster einer Zelle hat sich mit der Entwicklung von Hoch- 
durchsatz-Hybridisierungstechniken vollzogen. Bei diesen Ver- 
25 fahren wird auf einem sogenannten Microarray die Expressions- 
rate vieler 100 verschiedener Gene gleichzeitig getestet. Mit 
Hilfe dieser Verfahren ist es moglich, das Genexpressions- 
muster einer Zelle zu bestimmen. 



30 Dazu werden in der Regel die in der Zelle synthetisierten 

mRNA (messenger RNA) bestimmt. Die mRNA ist ein Zwischenpro- 
dukt bei der Translation des Gens zum Protein. Die mRNA ist 
somit eine Vorstufe bei der Bildung des Proteins und weist 
auf die Bildung des zugehorigen Proteins hin . Die zu untersu- 

35 chende Zelle wird zunachst isoliert. AnschlieJiend wird sie 
auf geschlossen . Durch geeignete Auf reinigungsschritte wird 
die mRNA aus der Zelle isoliert. Danach wird die mRNA mittels 
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der reversen Transkriptase in cDNA (complementary DNA) uber- 
setzt. Diese wird mit i. d. R. linearer PCR (polymerase chain 
reaction) amplif iziert . Die so gewonnene cDNA wird mit Hilfe 
von geeigneten Microarrays, z. B. DNA-Chips, qualitativ bzw. 
5 quantitativ analysiert. Mit modernen Microarrays konnen die 
Expressionsraten von 5.000 und mehr Genen gleichzeitig ver- 
messen werden. 

Aufgrund dieser verbesserten Techniken steht mittlerweile ein 
10 umfangreiches Wissen uber das menschliche Genom und Proteom 
sowie uber die Wechselwirkungen zwischen Proteinen und Genen 
bzw. Proteinen untereinander zur Verfugung. 

Von besonderer Bedeutung sind Datensatze, in denen die zeit- 
15 liche Abfolge von Genexpressionsmustern in einem Gewebe fest- 
gehalten ist. Als Beispiel sind sogenannte longitudinale 
Hybridisierungsstudien, d. h. zeitliche Abfolgen von Ge- 
nexpressionsmustern wahrend der Organdif f erenzierung im Rah- 
men der Embryonalentwicklung, zu nennen. Zeitauf geloste Ge- 
20 nexpressionsdaten existieren auch fur den Zellteilungs-Zyklus 
von Einzellern und sind auch fur komplexere Gewebe moglich. 

Neuronale Modellierung von Genom und Proteom (vgl . [5]) 

25 

Im Folgenden wird allgemein das Modellierungsprinzip geschil- 
dert. Grundzuge sind aus [5] bekannt. 

Das Grundprinzip besteht in der Herstellung einer Aquivalenz- 
30 beziehung zwischen dem f unktionellen Netzwerk des Genoms und 
Proteoms einerseits und dem Neuronalen Netz des menschlichen 
Gehirns andererseits, die beide stark vernetzte ruckgekoppel- 
te Systeme darstellen. 

35 Im Folgenden wird das Neuronale Netz des menschlichen Gehirns 
in einigen Grundziigen anhand von Fig. 2 geschildert. 
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Im menschlichen Gehirn befinden sich etwa hundert Milliarden 
Nervenzellen 20 (Neuronen) , die jede mit Zehntausenden ande- 
ren Nervenzellen 20 Inf ormationen austauscht. Die Informatio- 
nen gelangen von einem Neuron 20 liber das zu jedem Neuron ge- 
horende Axon 22 zu einem anderen Neuron 20. Jedes Neuron hat 
genau ein Axon, urn Inf ormationen an andere Neuronen zu sen- 
den. In seinem weiteren Verlauf verzweigt sich ein Axon typi- 
scherweise etwa eintausend Mai, so dass ein Neuron 20 uber 
sein Axon 22 Inf ormationen an etwa eintausend andere Neuronen 
20 senden kann. 

Zum Empfangen von Inf ormationen verfugen Neuronen 20 uber 
Dendriten 24. Das die Information tragende Axon 22 ist iiber 
eine Synapse 2 6 mit dem Dendriten 24 verbunden. fiber diese 
Synapse 2 6 gelangt die Information vom Axon 22 in den Dendrit 
24 und damit vom aussendenden Neuron 20 zum nachgeschalteten 
Neuron. An einem einzigen Dendriten 24 konnen Tausende bis 
Hunderttausende Axons 22 bzw. Synapsen 2 6 angreifen, so dass 
ein nachgeschaltetes Neuron 20 Signale von vielen 1000 vorge- 
schalteten Neuronen empfangen kann . 

Im Folgenden wird auf Fig. 3 Bezug genommen, die die Potenzi- 
ale innerhalb eines Dendriten bzw. eines Neurons als Funktion 
der Zeit darstellt. Der Inf ormationsaustausch zwischen den 
Neuronen 20 geschieht in Form von Aktionspotenzialen (Spikes) 
30, die jedes Neuron 20 iiber sein Axon 22 aussendet. Die Spi- 
kes evozieren in den nachgeschalteten Neuronen 20 erneut Sig- 
nale, die sogenannten postsynaptischen Potenziale (PSPs) 32. 
Die Groiie dieser PSPs hangt von der Ubertragungsstarke oder 
dem synaptischen Gewicht w der jeweiligen Synapse ab. 

Im Folgenden wird auf Fig. 4 Bezug genommen. Fig. 4 zeigt ei- 
nen Dendriten 24, an den eine erste Synapse 2 6 koppelt . An 
diese erste Synapse 26 ist eine zweite Synapse 36 gekoppelt. 
Diese zweite Synapse 36 wird als modulator ische Synapse be- 
zeichnet. Bezeichnen wir, in Anlehnung an Fig. 3, mit PSP das 
postsynaptische Potenzial 32, dass sich im Dendriten 24 auf- 
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grund der Einwirkung der ersten Synapse 2 6 in Abwesenheit der 
modulatorischen Synapse 36 ausbilden wtirde, so kann dies dar- 
gestellt werden durch 

5 PSP = w-e{t) , 

wobei w, wie oben definiert, das synaptische Gewicht der ers- 
ten Synapse 26 wiedergibt, und e(t) den zeitlichen Verlauf des 

postsynaptischen Potenzials 32 in einer geeigneten Normierung 

10 wiedergibt. 

Greift an die erste Synapse 2 6 zusatzlich die modulatorische 
Synapse 36 an, so ergibt sich im Dendriten 24 ein modifizier- 
tes postsynaptisches Potenzial PSP', das durch einen mul- 
15 tiplikativen Term act ausgedriickt werden kann: 

PSP' = act*we{t)=act-PSP . 

Dabei bezeichnet act die Aktivitat der modulatorischen Sy- 
20 napse 36. 

Beispielsweise haben im Zentralnervensystern, also im Gehirn 
und Ruckenmark, dopaminerge Synapsen modulatorischen Charak- 
ter . 

25 

Die neuronale Aktivitat jedes Neurons, d. h. die Anzahl der 
pro Zeiteinheit ausgesendeten Spikes, ergibt sich - verein- 
facht gesprochen - durch eine nichtlineare und zeitlich 
nichtlokale Funktion aller einlaufenden postsynaptischen Po- 
30 tenziale. Uberschreitet diese Funktion einen bestimmten 

Schwellenwert, so wird ein Spike 30 ausgelost und liber das 
Axon 22 ubertragen. 

Somit stellt das biologische Neuronale Netz des Gehirns ein 
35 komplexes, nichtlineares System dar, das obendrein eine hohe 
Vernetzungsdichte aufweist. Fur die modellhafte Beschreibung 
dieses Systems hat die Neuroinf ormatik in den vergangenen 
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Jahren leistungsf ahige Theorien und Algorithmen entwickelt 
( z . B . Kompartimentmodelle, Spike-Reponse-Modelle, Mean-Field- 
Modelle, Multimodulare Neurokognitive Modelle, Bayes-Belief- 
Netze) . 

5 

Diese Theorien bzw. Gleichungen entsprechen in ihrer Struktur 
den weiter oben fur die Reaktionskinetik abgeleiteten Glei- 
chungen. Daher kann das regulatorische Netzwerk aus Genom und 
Proteom einer Zelle wie folgt auf ein aquivalentes Neuronales 
10 Netz abgebildet werden : 

- Ein Gen A des Genoms (hier verstanden als diejenige Kombi- 
nation von Exons, die eindeutig ein Protein codiert) sowie 
das zugehorige Protein A werden mit einem Neuron A des aqui- 

15 valenten Neuronalen Netzes identif iziert . Da in Genexpressi- 
onsmustern nur die mRNA bzw. cDNA quantitativ analysiert 
wird, kann auf der Ebene der Genexpressionsmuster allein oh- 
nehin nicht zwischen Genen und Proteinen unterschieden wer- 
den. 

20 

- Die Expressionsrate eines Gens A wird als nichtnegative Ak- 
tivitat, z. B. die Spikerate des Neurons A auf gef asst. 

- Wirkt ein Protein B regulatorisch auf ein nachgeschaltetes 
25 Gen A, so enthalt das aquivalente Neuronale Netz eine synap- 

tische Verbindung von Neuron B zu dem aquivalenten nachge- 
schalteten Neuron A. 

- Die Art der regulatorischen Wirkung (verstarkend oder inhi- 
30 bitorisch) wird im Neuronalen Netz durch das Vorzeichen und 

die Starke des zugehorigen synaptischen Gewichts angegeben. 

- Eine posttranslationale Modifikation eines Proteins durch 
ein anderes Protein, in Fig. 1 z. B. die Modifikation von 

35 Protein B durch Protein D, entspricht der Wirkung einer modu- 
latorischen Synapse im Zentralnervensystem. Modulatorische 
Synapsen werden in kiinstlichen Neuronalen Netzen durch synap- 
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tische Verbindungen mit multiplikativer Wirkung auf andere 
Synapsen beschrieben. Die aquivalente Wiederspiegelung einer 
posttranslationalen Modifikation des Proteins B durch das 
Protein D ist also eine synaptische Verbindungen mit mul- 
5 tiplikativer Wirkung von Neuron D zu Neuron B. 

- Externe Signale werden mit Inputknoten des aquivalenten 
Neuronalen Netzes identif iziert . 

10 Zwischen genetischen und Neuronalen Netzen kann die beschrie- 
bene Aquivalenzbeziehung hergestellt werden. Das dynamische 
Wechselspiel von Genen und regulatorischen Proteinen wird da- 
mit durch ein dynamisches Neuronales Netz modelliert. 

15 Zu den geeigneten neuronalen Algorithmen zahlen Netze spiken- 
der Neuronen aber auch Mean-Field Modelle, die den expliziten 
Zeitablauf der Signaliibertragung zwischen den Neuronen durch 
die explizite Beschreibung der pos-tsynaptischen Potenziale 
berucksichtigen. Sie erlauben die Modellierung der zeitlichen 

20 Entwicklung der neuronalen Aktivitaten im Netz als Ergebnis 
externer Stimulation oder intrinsischer Aktivitat. 

Der zeitliche Verlauf der Konzentrationen, der sich durch die 
Reaktionskinetik zwischen den beteiligten Molekiilen (z. B. 
25 zwischen regulatorischem Protein und DNA-Promotor ) ergibt, 

wird also durch den Zeitverlauf der Aktivitaten der Neuronen 
ersetzt, so dass das resultierende Netzmodell zur Simulation 
der zeitlichen Entwicklung von Genexpressionsmustern herange- 
zogen werden kann. 

30 

Fur ein solches Neuronales Netz konnen die neuronalen Aktivi- 
taten iiber die Zeit erschlossen werden. Da die neuronalen Ak- 
tivitaten den Genexpressionsmustern entsprechen, konnen beide 
miteinander verglichen werden. Das Neuronale Netz entspricht 
35 einem simulierten Genexpressionsmuster . 



WO 2005/026722 



PCT/EP2004/051835 



15 

Ziel der Modellierung ist es, das dem Expressionsverlauf 
zugrurideliegende regulatorische Netzwerk zu ermitteln, also 
die Frage zu beantworten: "Welche neuronale Vernetzungsstruk- 
tur mit welchen Gewichten und Reaktionskonstanten ist konsi- 
5 stent mit dem beobachteten Genexpressionsverlauf?" 

Zur Beantwortung dieser Frage wird das Netz mit einem an 
Strukturlernen angelehnten Verfahren trainiert: Es wird ver- 
sucht, das beobachtete Verhalten mit so wenigen regulatori- 
10 schen Verbindungen wie moglich aber auch so gut wie moglich 
zu erklaren, also das einfachste mit den Daten konsistente 
Modell zu finden. 

Ein bevorzugtes Optimierungsverf ahren minimiert die totale 
Abweichung zwischen gemessenen und simulierten Genexpressi- 
onsmustern unter Verwendung eines "sparse priors ", also einer 
Zusatzbedingung, die die Koexistenz vieler Verbindungen mit 
kleinen Gewichten zugunsten weniger regulatorischer Verbin- 
dungen bestraft. Eine Moglichkeit, einen solchen "sparse pri- 
or" zu realisieren, ist der Fachwelt bekannt . 



15 



20 



Kreuzvalidierung und statistische Optimierung erlauben die 
Abschatzung der Eindeutigkeit der Losung sowie ihrer Vorher- 
sagekraf t (Generalisierungsf ahigkeit) . 

25 

Dem fertig trainierten Netz konnen anhand der Verbindungs- 
struktur der Neuronen kausale Zusammenhange zwischen Genen 
aber auch die Rolle verschiedener Gene entnommen werden. So 
weist ein asymmetrisches Gewicht nur von Gen B zu Gen A dar- 

30 auf hin, dass Gen B Gen A reguliert . Gleichzeitig konnen im 
Modell verschiedene Gene oder regulatorische Verbindungen 
kunstlich ab- oder angeschaltet und die Auswirkungen auf die 
Genexpressionsmuster mit dem Ziel quantif iziert werden , die 
Drsache (n) krankheitsbedingter Veranderungen der Genexpressi- 

35 on zu identif izieren (sogenannte inverse Modellierung) . 
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Im Rahmen dieses Dokuments sind folgende Verof f entlichungen 
zitiert : 

[1] Zell, A., "Simulation Neuronaler Netze", S. 35 bis 51, 
S. 55 bis 86, Addison-Wesley Longman Verlag GmbH, 1994, 
3. Unveranderter Nachdruck, R. Oldenbourg Verlag, ISBN 
3-486-24350-0, 2000 

[2] Zell, A., "Simulation Neuronaler Netze", S. 87 bis 114, 
Addison-Wesley Longman Verlag GmbH, 1994, 3. Unveran- 
derter Nachdruck, R. Oldenbourg Verlag, ISBN 3-486- 
24350-0, 2000 

[3] J. J. Binney, N. J. Dowrick, A. J* Fisher, M. E. J. New- 
man, "The Theory of Critical Phenomena", Kap. 6: Mean- 
Field Theory, Clarendon Press Oxford 1992 

[4] C. Koch, I. Segev (Hrsg) , "Methods of Neural Modeling: 
From Ions to Networks", Kap 13: D. Hansel and H. Sompo- 
linsky: "Modeling Feature Selectivity in Local Circu- 
its", MIT Press, Cambridge, 1998 

[5] Patentanmeldung PCT/DE02/03381 
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Patentanspruche 

1. Verfahren zum Identif izieren pharmazeutischer Targets 
mit folgenden Schritten: 

a) eine Mehrzahl von Genexpressionsmustern gleichartiger 
Zellen wird bestimmt, wobei jeweils die Expressionsrate der 
Gene der Zelle bestimmt wird; 

b) die Mehrzahl der Genexpressionsmuster wird derart be- 
stimmt, dass der zeitliche Verlauf der Genexpressionsmuster 
der Zelle zumindest teilweise rekonstruiert werden kann; 

c) ein dynamisches Modell des regulatorischen Netzwerks 
aus Genom und Proteom der Zelle wird gebildet, indem ein 
aquivalentes Neuronales Netz in der folgenden Weise gebildet 
wird: 

i) ein Gen des Genoms sowie das zugehorige Protein 
werden durch ein Neuron (20) des aquivalenten Neuronalen Net- 
zes reprasentiert; 

ii) die Expressionsrate eines Gens wird durch eine 
nichtnegative Aktivitat des aquivalenten Neurons reprasen- 
tiert; 

iii) die regulatorische Wirkung eines Proteins auf 
ein Gen wird durch eine synaptische Verbindung (26) von dern 
zum Protein aquivalenten Neuron zu dem zum Gen aquivalenten 
Neuron reprasentiert; 

iv) die Art einer regulatorischen Wirkung (verstar- 
kend oder inhibitorisch) wird im Neuronalen Netz durch das 
Vorzeichen und die Starke eines zugehorigen synaptischen Ge- 
wichts (w) reprasentiert; 

d) das aquivalente Neuronale Netz wird mit den bestimmten 
Genexpressionsmustern verglichen und an diese angepasst; 

e) aus dem angepassten Neuronalen Netz wird das regulato- 
rische Netzwerk der untersuchten Zelle erschlossen. 

2. Verfahren nach dem vorhergehenden Anspruch, 
dadurch gekennzeichnet, 

dass eine posttranslationale Modifikation eines ersten 
Proteins durch ein zweites Protein durch eine synaptische 
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Verbindungen (36) mit multiplikativer Wirkung vom zweiten 
Neuron zum ersten Neuron reprasentiert wird. 

3. Verfahren nach einem der vorhergehenden Anspriiche, 
dadurch gekennzeichnet, 

dass ein externer Einfluss auf das regulatorische Netz- 
werk durch einen Inputknoten des aquivalenten Neuronalen Net- 
zes reprasentiert wird. 

4. Verfahren nach einem der vorhergehenden Anspriiche, 
dadurch gekennzeichnet, 

dass bei der Anpassung des aquivalenten Neuronalen Netzes 
an die bestimmten Genexpressionsmuster solche Netze bevorzugt 
werden, die einen geringeren Grad an Vernetzung aufweisen. 

5. Anordnung zum Identif izieren pharmazeutischer Targets 
mit Mitteln zum Bestimmen einer Mehrzahl von Genexpressi- 

onsmustern einer Zelle, wobei jeweils die Expressionsrate der 
Gene der Zelle bestimmt wird, wobei die Mittel derart ausge- 
bildet sind, dass der zeitliche Verlauf der Genexpressions- 
muster der Zelle zumindest teilweise rekonstruiert werden 
kann; 

mit Mitteln zum Bilden eines dynamischen Modells des re- 
gulatorischen Netzwerks aus Genom und Proteom der Zelle, in- 
dent ein aquivalentes Neuronales Netz in der folgenden Weise 
gebildet wird: 

i) Ein Gen des Genoms sowie das zugehorige Protein 
werden durch ein Neuron (20) des aquivalenten Neuronalen Net- 
zes reprasentiert; 

ii) die Expressionsrate eines Gens wird durch eine 
nichtnegative Aktivitat des aquivalenten Neurons (2 0) repra- 
sentiert; 

iii) die regulatorische Wirkung eines Proteins auf 
ein Gen wird durch eine synaptische Verbindung <2 6) von dem 
zum Protein aquivalenten Neuron zu dem zum Gen aquivalenten 
Neuron reprasentiert ; 
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iv) die Art einer regulatorischen Wirkung (verstar- 
kend Oder inhibitorisch) wird im Neuronalen Netz durch das 
Vorzeichen und die Starke eines zugehorigen synaptischen Ge- 
wichts (w) reprasentiert; 
5 mit Mitteln zum Vergleichen des aquivalenten Neuronalen 

Netzes rait den bestimmten Genexpressionsmustern; 

mit Mitteln zum Anpassen des aquivalenten Neuronalen Net- 
zes an die bestimmten Genexpressionsmuster; und 

mit Mitteln zum Erschlieften des regulatorischen Netzwerks 
10 der untersuchten Zelle aus dem angepassten Neuronalen Netz. 
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